Dissipative Visco-plastic Deformation in Dynamic Fracture: 
Tip Blunting and Velocity Selection 
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Dynamic fracture in a wide class of materials reveals "fracture energy" T much larger than the 
expected nominal surface energy due to the formation of two fresh surfaces. Moreover, the fracture 
energy depends on the crack velocity, r = r(i;). We show that a simple dynamical theory of visco- 
plasticity coupled to asymptotic pure linear-elasticity provides a possible explanation to the above 
phenomena. The theory predicts tip blunting characterized by a dynamically determined crack tip 
radius of curvature. In addition, we demonstrate velocity selection for cracks in fixed-grip strip 
geometry accompanied by the identification of F and its velocity dependence. 



Introduction: The dynamics of rapid fracture are 
still posing a great challenge to theoretical physics [i|. 
A major difficulty in understanding the rich experimen- 
tal phenomenology is the coupling between very different 
scales; while fracture is ultimately driven by the release 
of the large scales elastic energy, the detailed deformation 
on the small scales near the moving crack edge determines 
both the stress conditions and energy dissipation that 
control the rate of growth and its direction. Since the 
near crack edge deformation is extremely complicated, it 
is common to lump all the near crack edge physics into a 
phenomenological unknown fracture energy, leaving the 
theory with a crucial missing ingredient. To better under- 
stand these issues, consider a semi-infinite crack forming 
in an infinite strip under fixed-grip boundary conditions. 
As the crack proceeds at velocity v, a constant energy per 
unit advance per unit time, J, is released from the elastic 
fields. In the classical theory of Linear Elasticity Frac- 
ture Mechanics one assumes that J — vT where F is the 
fracture energy, putatively a constant, equal to the en- 
ergy 27s of forming two fresh surfaces. In fact, in a large 
variety of materials, from steel to brittle plastics, the 
measured fracture energy F can be orders of magnitude 
larger than 27^. Moreover, in many instances the mea- 
sured fracture energy depends on the velocity, F = F(w). 
This is important; cracks tend to reach a stationary ve- 
locity V that grows as a function of the loading and this is 
possible only if F depends on v. It is tempting to assert 
that some kind of plastic processes may be responsible 
for this kind of findings. Straightforward plasticity the- 
ory assumes that the stress fields near the tip of a crack 
are tamed by plastic deformations and remain at a value 
close to the yield stress Sy of the material. Alas, this 
value of stress is assumed to be smaller than the level of 
stress necessary for overcoming the cohesive forces that 
bind material together. Without a proper theory of the 
dynamical processes occurring near the crack tip it re- 
mains very difficult to explain all these observations in a 
coherent and self-consistent manner. 

Theoretical framework: a promising formulation of 
visco-plasticity was proposed recently by Falk and Langer 
[J. As in all approaches to plasticity, one adds to the 
usual stress and elastic strain tensors aij and efj , a plas- 
tic component of the strain tensor, denoted as a trace- 



less symmetric tensor ef-. Essential to this theory are 
the highly localized regions called the "shear transfor- 
mation zones" (STZ). These dominate the plastic events 
and their presence is carried by two "internal state" fields 
U- One is a trace-less symmetric tensor field A that is 
related to the orientations of these zones and acts as a 
"back stress" ; the second is a scalar field A which stands 
for the density of these STZ. We employ here the tenso- 
rial quasi-linear version of this model, which is given by 
the following three coupled equations [J : 
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Here Sij is the 2D deviatoric stress tensor, Sij = aij — 
^TiaSij, r is a typical time scale and A is a parameter of 
dimension inverse stress that measures the plastic strain 
rate sensitivity to stress. This set of equations and fur- 
ther developments were studied in considerable detail in 
other contexts p, ^ . In this Letter we study the impli- 
cations of this model for dynamic fracture. The reader 
should note that in principle the time scale r may de- 
pend strongly on Sy , maybe even diverging for \sij \<$^Sy. 
With such divergence one can connect smoothly to elas- 
tic behavior far away from the crack tip. In our analysis 
below we take r to be constant, but on the other hand al- 
low the density A of STZ to decay to its pre-deformation 
small value far away from the crack tip. In this we differ 
from previous treatments Q where A was taken to be 
unity. We will show below that the dynamics of A are 
interesting and important, in addition to allowing us a 
smooth coupling to elasticity. 

Consider then a semi-infinite crack in a 2D infinite strip 
of a given width, loaded in fixed-grip boundary conditions 
such that the energy release rate G = J/v is experimen- 
tally controlled J3|- The crack advances steadily at a 
velocity v (to be found self-consistently) . The effect of 
the plastic response should be strongest near the tip of 
the crack, while far away from the tip linear elasticity 
can be safely employed. Like in |y] we expect the crack 



tip radius of curvature R to adjust itself dynamically, al- 
lowing the crack tip to blunt. The degree of blunting 
will be connected self-consistently to the steady velocity 
V '9'|. While we adopt below a number of approxima- 
tions, we insist on the conservation of energy in the sense 
that the stored elastic energy is mainly dissipated near 
the blunted tip (a feature that appears missing in |y]), 
though part of it can be locked in the crack wake. Sen- 
sitive to the experimental finding that r3>27s, we will 
assert that the dominant contribution to T comes from 
the plastic work which we assume to include both the 
energy dissipated and the residual energy locked in the 
wake [31; this assumption translates to 
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where A is the area in which the plastic dissipation is not 
vanishingly small. This equation should be considered 
together with the kinematic relation at the crack tip {9 — 
0), 
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where we have assumed that the plastic strain rate dom- 
inates the elastic one. 

In order to solve Eqs. Q-© we need an explicit so- 
lution for Sij. The solution for the deviatoric stress is 
not independent of the visco-plastic constitutive model 
in Eqs. C3l-®- The coupling is via the relation of the 
deviatoric stress to the elastic part of the strain tensor, 
which needs to add up to the plastic part such that the 
total strain satisfies the usual compatibility relation ^. 
We thus expect Sy to deviate from the usual linear-elastic 
behavior. To quantify this deviation we introduce the 
functions hij{R/r,6,v) (no summation convention) 



Before proceeding to our results, we recall a different 
line of investigation taken in JSJ. There the near sharp 
tip stress distribution of Eq. , with hij — 1 , was used 
to compute the plastic strain rate using a phenomenolog- 
ical visco-plastic constitutive law. Then the dissipation 
rate associated with the plastic strain rate was computed 
and added to an unknown velocity-independent dissipa- 
tion to be equated to the energy release rate G. The two 
approaches, ours and the one presented in [g , are approx- 
imate in nature and even share some of the underlying 
assumptions, but in fact are very different. Indeed, we 
use elastic stress solutions inside the plastic zone as in 
|3 and also demand explicitly energy conservation, but 
the similarity ends here. We consider the sharp tip as- 
sumption unjustified and by asserting that the crack tip 
advances mainly by plastic flow, arrive at Eq. jSJl that 
has no analog in [g; as a result the unrealistic stress 
singularity is cut off. Moreover, in contrast to j8| we 
assert that the plastic dissipation is dominant. Finally, 
we incorporate a dynamical rate- and- state formulation of 
visco-plasticity that is fundamentally different from the 
phenomenological constitutive law employed in 8] . 

Analysis: To analyze the various regimes of solutions 
we non-dimensionalize the equations by measuring stress 
in units of Sy, time in units of t and velocity in units of 
Cs- Having only one typical time scale, we expect that 
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hind the tip and Y^j 
The unknown functions h 



where r= a/(x — ut)^ + y^, 6 = ta.n~^ {y/{x — vt)) are po- 
lar coordinate system whose origin is at distance R be- 
are known universal functions Q- 
ij incorporate both a dynamic 
effect {sij appears now also in the time derivative of the 
compatibility relation) and a geometric effect (tip blunt- 
ing) on the stress. Clearly, these functions approach 
unity as R/r —^ and the remaining part is the sharp-tip 
linear-elastic asymptotic solution with the stress inten- 
sity factor Ki{v) (3 



In light of Eq. JSJ, sufficiently near the tip the dimen- 
sionless plastic strain rate is also of 0(1). In seeking 
solutions we keep this physical constraint in mind. 

• First regime of solutions, apparently not physical 

An analytically tractable parameter regime is Xsy <^ 1. 
Observe Eq. JT]), and note that the LHS is of 0(1). The 
hi,j{R/r, 0, v) , (6) fact that Xsy < 1 renders Sij/sy > 1, sending Aij to be 

■€i 1, according to Eq. jSJ). It also follows that A ^ 1 
very rapidly. In this limit the set of equations Q-® 
simplifies in the region of interest near the tip to 



A 

-Si 



(9) 



K,{v) 



EG 



A^{v/cs,v) 



(7) 



Here A^ is a known universal function [3, and Cg, E 
and V are the shear wave speed, Young's modulus and 
Poisson's ratio respectively. Below we approximate the 
functions hij in order to integrate Eqs. Q-® and find 
self-consistently v and R according to Eqs. ® and Q. 



Therefore the equations for the components of e^- de- 
couple in this limit, identifying A, as expected, with the 
visco-plastic analog of an inverse shear modulus. To pro- 
ceed, we assert that the predominant w-dependence in Eq. 
Q is carried by K-^{v). This is justified since E^ are qual- 
itatively the same for velocities below the Yoffe threshold 
[3 and the "u-dependence in hij can be neglected accord- 
ing to 7', 'q. We can substitute now Eq. ® in Eqs. H4I5|) 
and use the approximate Sij of (jHl we end up with 
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Note that this solution has interesting aspects: we find 
a self-consistent steady velocity that increases linearly 
upon increasing G. This implies that r(w)~t;. This solu- 
tion is realized self-consistently by having also R increas- 
ing with G. Since A/ — > 1 with zero slope when i; -^- 0^ 
[31, R starts linearly with G and later turns nonlinear 
according to the dependence of Aj on G. Unfortunately, 
this solution allows very large stress concentrations near 
the tip, i.e. Sij/sy ^ 1 and as we are unaware of experi- 
mental observations indicating such high stress levels, we 
propose that the limit Xsy <C 1 is not physically relevant. 

• Second regime of solutions, apparently physical 

The range of parameters that furnishes an interesting so- 
lution is Xsy ^ 1. In this range the constitutive equations 
(E3)-® are fully coupled and the "internal state" fields 
A and A play an important role in the dynamics. To see 
this we first note that in order to satisfy Eq. in the 
near tip region the combination XAsij should be of 0(1). 
According to Eqs. (^ and Q this condition can be re- 
alized with A < 1 and Sij ^ Sy. In that case, according 
to Eq. (O , the "back stress" Ay is comparable to XAsij 
and indeed resists plastic deformations. Thus for A.Sj,>l 
we expect a solution with Sij ^ Sy, where A does not 
saturate rapidly and the dynamics of A are important. 
To shed more light on the nature of this solution and to 
substantiate the first solution, we consider the equations 
numerically. 

Numerics: For the sake of numerical solutions we 
approximate the solution for Sij by allowing the visco- 
plastic processes to blunt the tip and then solve approxi- 
mately for the stress field up to the blunted tip via quasi- 
static linear elasticity. This approach inevitably violates 
compatibility, but is expected to yield reasonable approx- 
imations as shown in J6|. Thus the main effect of the 
visco-plastic processes on the stress is via the introduc- 
tion of the length scale R near the tip, hence neglecting 
all the w-dependence except for K^{v), as was explained 
above. In that case, outside the crack we solve for the 
stress field with boundary conditions 
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for a finite range of angles (say, —tt/3<9<tt/3). Since r 
cannot go to zero, we can add to the usual r~^/^ solution 
more "divergent" powers; explicitly, we add terms to (Ty 
of the form 
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The pre-factors of the new terms can be selected such as 
to best approximate Eq. H12|l in the sense that relative 
to the non- vanishing component agg the ratios Icr^r /creel 
and \(Jr$/(Jee\ are bounded below 1%. This solution is 
naturally matched to the usual linear-clastic asymptotic 
solution far from the tip, while capturing qualitatively 
the elimination of the common sharp tip singularity near 
the tip. 




FIG. 1: Solutions for Xsy = 1. Upper panel: The tip radius 
R normalized by the first data point as a function of the nor- 
malized velocity v/cs- Lower panel: The fracture energy F 
normalized by the first data point as a function of the nor- 
malized velocity v/cs- 



At this point we solve Eqs. Q-© along the line 6 — 0. 
The calculation of the integrals in Eq. ^ should be im- 
plemented with care. In contrast to classical theories of 
plasticity, here there is no sharp boundary between the 
visco-plastic and elastic solutions: the deviatoric stress 
Sij and therefore e^- decay algebraically. Thus we should 
choose a reasonable cutoff for which the integral in our 
model converges. For that purpose, we compute the in- 
tegrand of Q along the line 6* = and in particular 
determine the position x = W at which the plastic work 
rate s^e^ reaches 2% of its maximal value at the tip. 
Using this upper limit of integration the value of the in- 
tegral becomes insensitive to the long tail of the plastic 
rate of work. We then estimate the 2D integral as the ID 
integral times aW, with a being a dimcnsionless factor 
representing the full contribution from the loading and 
unloading regions (where in the latter the model is in- 
applicable 1^). The model equations exhibit a solution 
when a is chosen in the range 0.17 — 0.25. Finally, the 
ID integral can be computed over time rather than space 
(since all the functions depend on x—vt). This allows us 
to integrate Eqs. Q in time, capturing the essential his- 
tory dependence of e^- on the internal state fields A and 
A. At this point, Eqs. Q-© can be regarded as defin- 
ing two surfaces in an abstract v — R space. We search 
graphically for the zero of the intersection line of these 
surfaces as parametrized by G. 



Following this procedure we solved for R and v for 
the two ranges of Xsy considered above. First, we ana- 
lyzed the model with Xsy — 0.01, v — 0.3, E/sy — 50, 
a = 0.25, efj{t = 0) = A^ (i = 0) = 0, and A(i = 0) = 0.01 
and compared the numerical solution to the analytic 
counterpart. We verified that for similar small initial 
values of A the results are rather insensitive to the exact 
value. The numerical solution matches Eqs. H10|l - Hllf) 
extremely well. Nevertheless, the dimensionless combi- 
nation G/syR in the numerical solution equals « 3 x 10'^. 
Taking PMMA [l^ as a representative for the relevant 
class of materials, this number implies R « 5nm. This 
length scale is much smaller than expected for PMMA, 
thus supporting our proposition that the range Xsy <C 1 
is not physically relevant. On the other hand, in the 
parameter range Xsy ^ 1 we verified that both "inter- 
nal state" fields A and A play an important role in the 
dynamics and the dependence of both v and _R on G is 
non-linear. The dependence of R and F on u for Xsy — 1 
is shown in Fig. ^ Numerically we have found that to a 
very good approximation 



F ex vAi{v) , R ex v/Ai{v) . 
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Since Ai appears in the theory always in the combination 
Aj/E (cf. Eq. iQ), we can write F ex vSyAj{v)/E, 
R ex vE/syAi(v) where we non-dimensionalized E by 
Sy . Rearranging these equations and keeping dimensions 
right we end up with the implicit relations 
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where /„ and fn are two dimensionless nonlinear func- 
tions computed from the numerical solution, reflecting 
the importance of the internal state variables A and A 
in the solution. Note that the solution for R in H15(l 
differs from the one provided in [g even in the qua- 
sistatic limit, reflecting the difference between the two ap- 
proaches, mainly in the use here of Eq. I^J . In this equa- 
tion G appears twice, once on the RHS and once through 
the stress intensity factor on the LHS, thus changing the 
way that E appears in the final solution. 

Finally, the dimensionless combination G/syR for this 
class of solutions is typically of 0(1). Applying the anal- 



ysis to PMMA we estimate R w 15/im, which is an ac- 
ceptable estimate for this class of materials. We propose 
therefore that the solutions in the range Xsy > 1 are phys- 
ically relevant for materials in which the main dissipation 
mechanism is visco-plasticity The stresses transmitted 
to the crack tip itself are moderately larger than Sy, sup- 
posedly sufficient to overcome the cohesive forces binding 
the material together, and the plastic zones are relatively 
thin and changing smoothly towards the linear elastic 
asymptotic fields. 

Summary: We have employed a theoretical frame- 
work of visco-plasticity to propose a direct measure of 
the fracture energy in dynamic fracture for a class of ma- 
terials where the fracture energy exceeds significantly the 
bare surface energy. The fracture energy is a function of 
the crack velocity, rationalizing the experimental obser- 
vations of velocity selection for cracks in a fixed-grip strip 
geometry. This selection is accompanied by the blunting 
of the tip, with a radius that is self-consistent with the 
velocity. The theory, within the stated approximations, 
predicts that both F and R are monotonically increasing 
functions of the velocity v in the considered range. In the 
physical range of parameters (Xsy > 1) the stress at the 
tip exceeds the yield-stress Sy, affording rupture of the 
cohesive bonds. In this range of parameters the plastic 
zone is relatively thin with respect to the tip radius. 

The STZ model as used here cannot be considered a 
final theory of visco-plasticity |5(. A rigorous coupling 
to elasticity and a consistent description of the wake of 
the crack are still missing here. Nevertheless we believe 
that the encouraging approximate results presented here 
are generic to any reasonable model that addresses the 
complex state of deformation near a moving crack tip. In 
future work one should invest effort in improved models 
where the approach can be tested in some quantitative 
detail. In such future work the questions of crack tip 
instabilities that were not addressed here at all should 
be directly addressed and explained. 
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